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1.  INTRODUCTION 


Tuberous  sclerosis  complex  (TSC)  is  an  inherited  disorder  characterized  by  the  appearance 
of  benign  tumors  in  the  brain,  heart,  lung,  kidney,  eyes,  and/or  skin.  The  disease  affects  as  many 
as  40,000  individuals  in  the  United  States  and  is  caused  by  mutations  in  one  of  two  genes,  TSCl 
(encoding  a  protein  called  hamartin),  and  TSC2  (encoding  a  protein  called  tuberin).  While 
research  of  the  past  two  decades  has  provided  a  wealth  of  information  regarding  the  genetic  origins 
of  TSC,  the  complexity  of  the  molecular  pathways  involved,  together  with  their  intricate  temporal 
and  spatial  interactions,  has  made  it  difficult  to  elucidate  the  full  molecular  pathology  of  the 
disease.  Fortunately,  recent  technical  and  conceptual  advances  with  respect  to  the  importance  of 
“genetic  buffering”  on  phenotypic  variation  have  provided  novel  avenues  of  exploration  regarding 
this  goal. 

Using  the  genetically  tractable  fission  yeast  as  a  model,  our  USAMRAA  funded  research 
has  sought  to  exploit  recent  advances  in  genetic  interaction  network  biology  to  analyze  previously 
uncharacterized  genes  that  modulate  the  phenotypic  effects  of  hamartin  and  tuberin  loss  of 
function  mutations.  Of  particular  interest  are  genes  displaying  negative  genetic  interactions  with 
either  tscl  or  tsc2.  Since  tumor  formation  in  TSC  patients  arises  from  loss  of  heterozygosity,  this 
characteristic  identifies  the  human  orthologs  of  these  genes  as  potential  therapeutic  targets  i.e. 
drugs  inhibiting  a  negative  interactor  would  presumably  suppress  only  the  growth  of  tumor  cells 
(which  bear  two  mutant  copies  of  the  affected  TSC  gene;  the  inherited  mutant  germline  copy,  and 
the  copy  affected  by  the  “second-hit”)  while  leaving  phenotypically  normal  cells  (carrying  only 
the  mutant  germline  copy)  unaffected. 

Our  research  of  the  past  year  has  clearly  identified  two  genes,  fft2  (encoding  a 
SMARCADl  family  ATP-dependent  DNA  helicase)  andypai  (encoding  a  PTPA  family  protein 
phosphatase)  as  excellent  candidates  for  continued  analysis.  While  deletion  of  either  gene  has 
little  phenotypic  effect  in  normal  cells,  their  loss  in  either  tscl  or  tsc2  mutant  backgrounds 
profoundly  inhibits  growth.  Thus,  inhibition  of  either  fft3  ox  ypal  may  represent  an  “Achilles’ 
heel”  of  cells  defective  in  hamartin  or  tuberin  function.  Importantly,  both  of  these  genes  have 
clear  orthologs  in  humans. 

In  addition  to  our  bench-top  analyses,  we  have  also  developed  a  formal  mathematical 
methodology  for  the  quantitation  of  genetic  buffering  strength  (this  was  necessitated  by  the  need 
for  a  system  that  could  be  used  to  unambiguously  compare  genetic  interaction  data).  This 
methodology  provides  a  simple,  general,  and  rigorous  mathematical  paradigm  with  which  to 
exactly  quantitate  the  buffering  strength  of  any  genetic  determinant.  The  methodology  was 
recently  published  in  the  peer-reviewed  “77ze  Mathematica  Journar  and  is  freely  available  to  the 
public  at  “www.mathematica-journal.eom/2015/03/the-quantitation-of-non-classical-buffering/”. 
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3.  ACCOMPLISHMENTS 


Our  research  accomplishments  of  the  last  year  are  described  below  in  the  context  of  the 
Specific  Aims  outlined  in  the  approved  Statement  of  Work  (SOW). 

Specific  Aim  1:  Establish  an  5.  pombe  model  of  TSC  by  replacing  the  endogenous  5.  pombe 
hamartin  and  tuberin  genes  with  human  TSCl  and  TSC2,  respectively. 

Task  1:  Create  tscl::ura4^  and  tsc2::ura4^  "base”  strains. 

Predicted  Completion  Date:  August  3  P',  20 1 4 

Progress:  Gene  deletion  cassettes  for  both  the  tscl  and  tsc2  genes  were  created  by  PCR 
amplifying  the  ura4^  selectable  marker  (together  with  the  loxP  and  loxM3  recombination 
sites)  from  the  pAWl  vector  using  primers  oTscl  1  (5’-tta  tea  atg  ctg  cca  aga  ctt  get  ate  agt 
ata  atg  teg  cat  agt  tgt  ata  tea  acg  ttg  act  ttg  cca  act  ttg  tac  gac  gga  tec  ccg  ggt  taa  tta  a-3’), 
and  oTscl2  (5’-aat  tat  ttt  ata  tgg  aat  gag  caa  gta  tgt  ttt  ate  ata  att  gac  cag  ttc  att  tea  agg  ace 
ttc  aaa  aat  ata  cct  acg  aat  teg  age  teg  ttt  aaa  c-3’),  or  oTscl3  (5’-tta  aga  gtt  cag  att  tgc  ttt  atg 
tgg  tta  ttc  tgc  tga  agg  tec  taa  ttt  att  gac  gtt  gaa  aaa  taa  agg  cca  cat  age  gga  tec  ccg  ggt  taa  tta 
a-3’),  and  oTscl4  (5 ’-ata  aaa  aaa  att  aat  taa  tga  tgg  caa  ggc  aca  ate  gta  ate  aat  ctt  tta  att  tag 
gac  ttt  tta  tat  gee  ctt  atg  gcg  aat  teg  age  teg  ttt  aaa  c-3’),  respectively.  Strain  ED666 
{ura4-D18  leul-32  ade6-210  K^)  was  then  transformed  with  either  the  tscl  specific  or  tsc2 
specific  cassette.  Ura"^  transformants  were  isolated  and  the  respective  gene  deletions 
confirmed  by  colony  PCR  using  primers  oTscl5  (5 ’-atg  tgg  cag  act  acg  eta  tec  t-3’)  and 
oTsclV  (5’-atg  ctt  ccc  eta  att  cat  age  a-3’)  for  the  tscl  deletion,  and  oTscl6  (5’-agc  aac  eta 
cga  gag  gaa  gat  g-3’)  and  oTsclS  (5 ’-gcg  cat  aac  cct  ttc  tac  att  c-3’)  for  the  tsc2  deletion. 

Status:  Complete. 


Task  2:  Obtain  full  length  cDNA  clones  of  human  TSCl  (Accession  BC167824)  and 
TSC2  (Accession  BC 150300)  from  ThermoScientific,  PCR  amplification. 

Predicted  Completion  Date:  August  3  P‘,  20 1 4 

Progress:  Recent  advances  in  DNA  synthesis  technology  altered  our  strategy  with  respect 
to  the  original  SOW  document  (see  Section  5:  Problems/Changes).  Instead  of  purchasing 
cDNA  clones  and  using  standard  amplification/cloning  techniques,  it  proved  more 
economical  (in  terms  of  both  time  and  money)  to  synthesize  the  desired  sequences.  Full 
length  clones  were  thus  obtained  through  Genscript’s  DNA  synthesis  service.  All 
constructs  were  confirmed  by  DNA  sequencing. 

Status:  Complete. 
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Task  3:  Molecular  cloning  of  TSCl  and  TSC2  into  "exchange"  plasmids. 

Predicted  Completion  Date:  September  30*,  2014 

Progress:  As  part  of  the  service  provided  by  Genscript,  it  was  possible  to  have  the  synthetic 
sequences  cloned  directly  into  the  pAWSX  “exchange”  plasmids  upon  synthesis. 

Status:  Complete. 


Task  4:  Transform  exchange  plasmids  into  5.  pombe  strain  ED666,  Select  for 
clones  that  have  exchanged  the  ura4*  cassette  with  TSC1ITSC2.  Verify  genotypes  hy 
colony  PCR. 

Predicted  Completion  Date:  October  3P\  2014 

Progress:  The  tscl  ::ura4/tsc2:ura4  base  strains  were  transformed  with  the 
pAW8X-TSCl/pAW8X-TSC2  vectors  containing  the  human  TSCl  and  TSC2  genes.  As 
controls,  the  base  strains  were  also  transformed  with  the  pAW8X-tscl  and  pAW8X-tsc2 
vectors  containing  the  fission  yeast  tscl  and  tsc2  genes.  Cells  in  which  the  um4^  gene 
was  exchanged  with  TSCl/TSC2/tscl/tsc2  were  selected  by  growth  on  media  containing 
5-fiuoroorotic  acid  (a  drug  counter  selectable  to  Ura^  cells).  This  created  strains  in  which 
human  TSCl  or  TSC2,  or  fission  yeast  tscl  or  tsc2,  were  expressed  from  the  endogenous 
fission  yeast  tscl/tsc2  promoters  at  the  native  tscl/tsc2  loci.  These  strains  have  been 
denoted  as  tscl::TSCl^\  tsc2::TSC2^\  tscl ::tscl^P,  and  tsc2::tsc2^P.  Strains  bearing  both 
substitutions  were  then  created  by  crossing  tscl::TSCl^^  and  tsc2::TSC2^^  strains  of 
opposite  mating  type  and  screening  progeny  by  colony  PCR  to  identify  tscl::TSCl^^ 
tsc2::TSC2^^  cells.  A  complete  catalog  of  strains  created  during  the  course  of  this  work 
appears  in  Section  6.5.1,  Research  Material. 

Status:  Complete. 


Task  5:  Complementation  Assays, 

Predicted  Completion  Date:  December  3  P*,  2014 

Progress:  The  ability  of  human  TSC1/TSC2  to  complement  the  loss  of  the  fission  yeast 
tscl/tsc2  genes  was  determined  by  assaying  the  growth  of  wild-type,  tscl  A,  tsc2A, 
tscl::TSCl^\  tsc2::TSC2^\  tscl::TSCl^^  tsc2::TSC2^\  tscl::tscl^P,  and  tsc2::tsc2^P 
strains  on  minimal  media  containing  60  pg/ml  of  canavanine.  As  expected  tsclA  and 
tsc2A  base  strains,  unlike  wild-type  controls,  were  resistant  to  the  drug  and  were  able  to 
form  colonies.  Also  as  expected,  tscl::tscl^P  and  tsc2::tsc2^P  strains  behaved  as  wild  type 
and  were  unable  to  form  colonies  (indicating  that  the  recombinase  mediated  cassette 
exchange  methodology  correctly  integrated  the  synthetic  constructs).  Unfortunately, 
tscl ::TSC1^\  tsc2: :TSC2^\  as  well  as  tscl::TSCl^Usc2::TSC2^^  strains  were  able  to  form 
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colonies  (indicating  that  the  expression  of  the  human  orthologs  was  incapable  of 
complementing  the  loss  of  the  fission  yeast  tscl  and  tsc2  genes). 

Status:  Complete. 


Specific  Aim  2:  Create  "TSC  mutant  allele"  array. 

Task  1:  Construction  of  TSC  mutant  alleles  via  direct  synthesis  or  site-directed 
mutagenesis  techniques. 

Predicted  Completion  Date:  February  28*,  2015 

Progress:  As  discussed  above,  due  to  the  lack  of  complementation,  the  remainder  of  the 
project  will  be  conducted  using  the  fission  yeast  tscl  and  tsc2  genes.  To  this  end, 
site-directed  mutagenesis  was  employed  to  incorporate  a  total  of  five  mutations  into  the 
tsc2  gene.  Each  of  these  mutations  (G296E,  R927W,  N1 199S,  P1223E,  R1296P)  affects 
residues  analogous  to  mutations  identified  clinically.  New  England  Biolab’s  Q5 
Site -Directed  Mutagenesis  Kit  was  used  to  generate  pUC57-tsc2_G296E, 
pUC57-tsc2_R927W,  pUC57_Nl  199S,  pUC57_P1223E,  and  pUC57_R1296P  plasmids. 
DNA  sequence  analysis  confirmed  the  incorporation  of  the  desired  changes.  The  same 
strategy  will  be  used  to  create  a  similar  panel  of  tscl  mutant  alleles. 

Status:  In  Progress  (50%  Complete). 


Task  2:  Molecular  cloning  of  mutant  alleles  into  exchange  plasmids. 

Predicted  Completion  Date:  March  3  E',  2015 

Progress:  As  of  June  30*,  2015,  two  of  the  clones  (pUC57-tsc2_G296E,  and 
pUC57_Nl  199S)  have  been  cloned  into  the  pAW8X  exchange  plasmids  using  the  Xhol 
and  Sad  restriction  sites  and  standard  molecular  techniques.  Cloning  of  the  remaining 
constructs  is  in  progress. 

Status:  In  Progress  (20%  Complete). 


Task  3:  Transform  exchange  plasmids  into  S.  pomhe  strain  ED666,  Select  for 
clones  that  have  exchanged  the  ura4*  cassette  with  TSC1/TSC2,  Verify  genotypes  hy 
colony  PCR, 

Predicted  Completion  Date:  April  3  E',  2015 

Progress:  Upon  completion  of  Task  2  (above),  recombinase  mediated  cassette  exchange 
will  be  used  to  integrate  the  respective  mutant  tscl  and  tsc2  alleles  into  the  appropriate  base 
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strains. 


Status:  Pending  Completion  of  Task  2. 

Task  4:  Confirm  lack  of  complementation. 

Predicted  Completion  Date:  June  SO*'’,  2015 

Progress:  Upon  completion  of  Tasks  2  and  3  (above),  complementation  assays  based  on 
canavanine  resistance  will  be  performed  as  described  in  Specific  Aim  1,  Task  5. 

Status:  Pending  Completion  of  Tasks  2  and  3. 


Specific  Aim  3:  Intercross  strains  of  the  TSC  mutant  allele  panel  with  strains  of  the 
"putative  interactor"  panel. 

Task  1:  Intercross  "TSC  mutant  allele"  array  with  "interactor"  panel. 

Predicted  Completion  Date:  October  3  U',  2015 

Progress:  Control  experiments  in  which  the  fission  yeast  interactor  panel  was  crossed  to 
either  the  tsclA  or  tsc2A  gene  deletion  mutants  have  been  completed.  These  control 
experiments  are  crucial  in  that  they  will  provide  a  proper  baseline  of  comparison  for  future 
assays  conducted  with  the  TSC  mutant  allele  panels.  Due  to  the  confounding  effect  of 
auxotrophic  markers  present  in  the  interactor  panel  (derived  from  the  Bioneer  gene 
deletion  set),  it  was  necessary  to  construct  an  alternative  panel  of  strains  devoid  of 
secondary  genetic  alterations  (see  Section  5:  Problems/Changes).  Intercrossing  of  the 
tscl  and  tsc2  mutant  allele  panels  will  be  performed  upon  completion  of  Specific  Aim  2. 

Status:  In  Progress  (30%  Complete). 


Task  2:  Isolation  of  double  mutant  progeny.  Growth  Assays. 

Predicted  Completion  Date:  October  3  U',  2015 

Progress:  Control  assays  (in  triplicate)  comparing  the  growth  of  wild  type,  tsclA,  tsc2A, 
the  interactor  panel  mutants,  as  well  as  the  respective  double  mutants  have  been  completed. 
Briefly,  a  micromanipulator  was  used  to  array  single  cells  representing  each  genotype  (in 
quadruplicate)  upon  YES  agar  plates.  Plates  were  incubated  at  30°C  for  five  days.  Digital 
images  were  taken  at  24  hour  intervals.  These  experiments  identified  jft3  (a  SMARCADl 
family  ATP-dependent  DNA  helicase)  and  ypal  (a  PTPA  family  protein  phosphatase)  as 
strong  negative  interactors  of  both  the  tsclA  and  tsc2A  gene  deletions  (see  Section  4.1, 
Impact).  Representative  images  of  plates  at  day  4  for  one  biological  replicate  are  shown 
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below  (10-27  denotes  ihs  ypal A  mutant,  and  30-77  denotes  the  fftSA  mutant).  Growth 
assays  using  the  tscl  and  tsc2  mutant  allele  panels  will  be  performed  upon  eompletion  of 
Speeific  Aim  2. 


Status:  In  Progress  (30%  Complete). 
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Task  3:  Image  analysis  of  colony  size. 

Predicted  Completion  Date:  December  3  P',  2015 

Progress:  Image!  (imagej.nih.gov/ij/)  analysis  tools  were  used  to  determine  the  colony  size 
of  wild  type,  tsclA,  tsc2A,  each  interactor  panel  mutant,  and  the  respective  double  mutants 
(in  quadruplicate)  for  each  of  three  independent  biological  replicates.  Representative  data 
describing  the  interaction  between  tscl  and  either  ypai  or  fft3  is  shown  below. 


Status:  In  Progress  (30%  Complete). 

Specific  Aim  4:  Physiological  analysis  of  TSC  modifiers 

Task  1:  Test  effect  of  pharmacological  inhibitors  of  "negative"  interactors  (if 
available)  in  strains  expressing  wild-type  and  mutant  TSC  alleles. 

Predicted  Completion  Date:  February  28*,  2016 

Progress:  These  assays  will  be  performed  upon  completion  of  Specific  Aims  2  and  3. 
Status:  Pending  Completion  of  Specific  Aims  2  and  3. 
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Task  2:  Molecular/physiological/genetic  analysis  of  interactors.  Specific  analyses 
will  be  determined  based  on  the  molecular  identity  of  the  genes  identified  from  the 
interactor  panel. 

Predicted  Completion  Date:  June  SO*'’,  2016 

Progress:  These  assays  will  be  performed  upon  completion  of  Specific  Aims  2  and  3. 
Status:  Pending  Completion  of  Specific  Aims  2  and  3. 


4,  IMPACT 

4.1  Impact  on  the  Development  of  the  Principal  Disciplines  of  the  Project 

Our  ongoing  study  has  clearly  identified  jft3  (a  SMARCADl  family 
ATP-dependent  DNA  helicase)  and  ypal  (a  PTPA  family  protein  phosphatase)  as  strong 
negative  interactors  of  both  the  tsclA  and  tsc2A  gene  deletions.  While  deletion  of  either 
gene  has  little  phenotypic  effect  in  normal  cells,  their  loss  in  either  tscl  or  tsc2  mutants 
profoundly  inhibits  growth.  Thus,  the  chemical  inhibition  of  the  encoded  Fft3  or  Ypal 
proteins  (both  of  which  have  clear  orthologs  in  humans)  may  represent  a  novel  means  with 
which  to  specifically  inhibit  the  growth  of  cells  defective  in  hamartin  or  tuberin  function. 
In  other  words,  the  inhibition  of  either  Fft3  or  Ypal  may  represent  an  “Achilles’  heel”  of 
cells  diminished  in  tscl  or  tsc2  activity.  It  is  important  to  note  that  neither  of  these  genes 
have  been  identified  using  traditional  genetic  analysis.  Thus,  their  continued  study 
represents  a  novel  avenue  of  exploration  with  regards  to  the  molecular  pathology  of  TSC. 

4.2  Impact  on  other  Disciplines 

During  the  course  of  analyzing  our  genetic  interaction  data,  it  became  apparent  that 
there  was  no  established  method  with  which  to  formally  and  unambiguously  quantitate  the 
strength  of  “genetic  buffering”.  To  address  this  issue  we  developed  a  formal,  axiomatic 
methodology  -  using  the  Wolfram  Programming  Language  -  to  achieve  this  goal.  The 
method  is  described  in  our  recent  publication  “The  Quantitation  of  Non-Classical 
Buffering:  Applying  the  Formal  and  General  Approach  to  Problems  in  the  Biological 
Sciences  (included  in  this  report  as  Appendix  A  and  freely  available  at 
“www.mathematica-journal.eom/2015/03/the-quantitation-of-non-classical-buffering/”. 
While  we  will  use  this  methodology  to  quantitate  and  compare  the  genetic  interactions 
observed  between  the  interactor  panel  and  TSC  mutant  allele  array,  the  methodology  can 
also  be  applied  to  any  scenario  (and  in  any  discipline)  where  a  quantity  partitions  between 
two  compartments  or  states. 

4.3  Impact  on  Technology  Transfer 

Nothing  to  Report 

4.4  Impact  on  Society  Beyond  Science  and  Technology 

Nothing  to  Report. 
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5.  CHANGES/PROBLEMS 


5,1  Changes  in  Approach 

5.1.1,  Cloning  of  human  TSCl  and  TSC2 

Instead  of  purchasing  cDNA  clones  of  TSCl  and  TSC2  and  subsequently 
employing  standard  amplification  and  cloning  techniques  to  create  the  exchange 
plasmids,  we  chose  to  use  a  commercial  entity  to  synthesize  the  desired  sequences. 
Recent  advances  in  DNA  synthesis  technology  have  made  such  direct  DNA  synthesis 
approaches  the  most  economical  choice.  Full  length  clones  were  thus  obtained  through 
Genscript’s  DNA  synthesis  service. 

5.1.2,  TSC  mutant  allele  panel 

Unfortunately,  human  TSCl  or  TSC2  were  unable  to  complement  the  loss  of 
their  fission  yeast  counterparts  when  expressed  at  the  native  tscl  or  tsc2  loci.  While 
disappointing,  this  was  an  eventuality  we  previously  considered.  As  discussed  in  the 
original  project  narrative  (page  4,  paragraph  1),  the  project  will  continue  using  the 
fission  yeast  hamartin/tuberin  genes  instead  of  their  human  orthologs.  This  is  to  say 
alleles  of  the  fission  yeast  tscl  and  tsc2  genes  -  bearing  mutations  affecting  residues 
analogous  to  those  identified  clinically  -  will  be  employed. 


5,2  Actual  or  Anticipated  Problems 

5,2,1  Confounding  Effect  of  Auxotrophic  Markers 

Due  to  the  confounding  effect  of  auxotrophic  markers  present  in  the  interactor 
panel  (derived  from  the  Bioneer  gene  deletion  set),  it  was  necessary  to  construct  an 
alternative  panel  devoid  of  secondary  genetic  alterations.  This  delayed  the  initiation  of 
Specific  Aims  2  and  3.  While  tedious  and  time  consuming  an  alternative  interactor 
panel,  devoid  of  any  background  auxotrophic  markers,  has  now  been  constructed  and 
will  be  employed  during  the  remainder  of  the  study. 


5,2,2  Personnel 

The  initial  proposal  identified  graduate  student,  Bidhan  Chakraborty,  as 
associated  with  the  project.  Unfortunately,  Mr.  Chakraborty,  upon  successfully 
defending  his  MSc  thesis,  returned  to  Bangladesh  to  be  married  and  ultimately  decided 
to  pursue  other  interests.  While  Mr.  Chakraborty  has  now  been  replaced  by  Mr.  Ryan 
Chevalier,  his  initial  departure  caused  minor  delays  with  respect  to  the  completion  of 
Specific  Aim  2.  These  manpower  issues  have  now  been  resolved  and  we  expect  no 
further  difficulties. 
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6.  PRODUCTS 


6,1  Publications,  Conference  Papers,  and  Presentations 

6.1.1  Journal  Publications 

Karagiannis  J.  The  Quantitation  of  Non-Classical  Buffering:  Applying  the 
Formal  and  General  Approach  to  Problems  in  the  Biological  Sciences.  The 
Mathematica  Journal  17,  2015.  (Published) 

6.1.2  Books  or  Other  Non-Periodical,  One-time  Publications 
Nothing  to  Report 

6.1.3  Other  Publications,  Conference  Papers,  and  Presentations 
Nothing  to  Report 


6,2  Websites  or  Other  Internet  Sites 

Nothing  to  Report 


6,3  Technologies  or  Techniques 

Nothing  to  Report 


6.4  Inventions,  Patent  Applications,  and/or  Licenses 

Nothing  to  Report 

6.5  Other  Products 

6,5.1  Research  Material  (Fission  Yeast  Strains) 

Genotypes  of  the  relevant  fission  yeast  strains  generated  during  the  course  of  this 
study  are  listed  below: 

•  tsclA  “base”  strain  (gene  deletion  with  flanking  loxP  and  loxM3  sites) 

•  tsc2A  “base”  strain  (gene  deletion  with  flanking  loxP  and  loxM3  sites) 

•  tscl::tscl^P 

•  tsc2::tsc2^P 

•  tscl::TSCl^^ 

•  tsc2::TSC2^^ 

•  tscl::TSCl^Jsc2::TSC2^^ 


10 


7.  PARTICIPANTS 


Name; 

Nearest  Person  Month  Worked: 
Contribution  to  Projeet: 


Jim  Karagiannis 
12 

Prineipal  Investigator 


Name;  Ryan  Chevalier 

Nearest  Person  Month  Worked  4 

Contribution  to  Projeet  Mr.  Chevalier  was  involved  in  the  routine  growth 

and  maintenanee  of  fission  yeast  cultures, 
molecular  cloning,  and  was  also  responsible  for  all 
image  analysis  of  genetic  interaction  data. 


8.  SPECIAL  REPORTING  REQUIREMENTS 

Not  Applicable 
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Appendix  A 


The  Mathematica®  Journal 


The  Quantitation  of  Non- 
classical  Buffering 

Applying  the  Formal  and  General  Approach 
to  Problems  in  the  Biological  Sciences 

Jim  Karagiannis 


A  formal,  axiomatic  conceptualization  of  buffering  action — 
generally  applicable  to  any  physical,  chemical,  or  biological 
process — was  first  presented  by  B.  M.  Schmitt  in  2005  [1 , 2], 
This  article  provides  a  series  of  tools  designed  to  aid  in  the 
application  of  these  concepts  to  both  classical  and  non-classical 
buffering  phenomena.  To  illustrate  the  utility  of  the  approach  in 
the  biological  sciences,  an  abstract  measure  of  the  magnitude  of 
“genetic”  buffering  associated  with  an  allele  of  the  gene 
encoding  the  heat  shock  protein  Hsp90  is  determined. 


■  Introduction 

Buffering  is  observed  when  a  parameter  changes  less  than  expected  in  response  to  a  given 
disturbance.  For  example,  if  a  strong  acid  is  added  to  an  aqueous  solution,  and  not  all  of 
the  added  ions  remain  free  (unbound),  then  the  solution  is  said  to  act  as  a  buffer. 
While  such  phenomena  have  been  observed  in  the  physical,  chemical,  and  biological  sci¬ 
ences  for  centuries,  a  formal  and  general  approach  for  their  quantitation  across  distinct  dis¬ 
ciplines  was  not  presented  until  recently  [1,2]. 

In  [1,  2],  B.  M.  Schmitt  presents  a  comprehensive  mathematical  framework  for  evaluating 
buffering  action — a  framework  that  is  applicable  to  any  scenario  in  which  a  quantity  parti¬ 
tions  between  two  compartments  or  states.  While  suitable  for  the  analysis  of  “classical” 
buffering  phenomena  (e.g.  acid-base  chemistry),  the  formalism  also  provides  a  simple 
means  with  which  to  quantitate  and  characterize  phenomena  that,  at  first  glance,  seem  to 
be  outside  the  buffering  paradigm — at  least  when  the  term  is  used  in  its  traditional  sense 
with  respect  to  the  homeostasis  of  physiological  parameters  (see  [2]  for  several  in-depth 
examples  of  the  formal  and  general  approach  applied  to  phenomena  spanning  diverse 
disciplines). 

This  article  provides  a  series  of  mathematical  tools  designed  to  facilitate  the  application  of 
Schmitt’s  paradigm  to  both  “classical”  and  “non-classical”  buffering  phenomena.  By 
inputting  the  transfer  function  t(x)  and  the  buffering  function  yS(x),  the  provided  code: 
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(1)  defines  the  four  buffering  parameters  t,  b,  T,  and  5;  (2)  provides  a  series  of  graphical 
outputs  describing  the  system;  and  (3)  calculates  the  buffering  angle  a,  thereby  classifying 
the  system  as  being  moderating/amplifying/inverting/non-inverting.  Thus,  the  burden  of 
computation  associated  with  the  approach  is  seamlessly  transferred  from  the  user  to  the 
Mathematica  kernel. 

Lastly,  in  order  to  illustrate  the  utility  of  the  approach  outside  traditional  disciplines,  I  use 
these  same  tools  to  calculate  the  magnitude  of  “genetic”  buffering  associated  with  an  al¬ 
lele  of  the  gene  encoding  the  heat  shock  protein  Hsp90.  In  this  way,  the  capacity  of  an  al¬ 
lele  to  buffer  phenotypic  variation  is  formally  and  quantitatively  determined  for  the  first 
time. 


■  An  Intuitive,  Linear  Example  Based  on  the  Partitioning  of 
Fluids  in  Communicating  Vessels 

To  illustrate  Schmitt’s  “formal  and  general  approach,”  consider  a  system  comprised  of  a 
pair  of  communicating  vessels  (i.e.  cylindrical  flasks  connected  by  a  small  tube).  As  de¬ 
scribed  in  [1],  the  total  volume  of  liquid  in  the  vessels,  x,  is  the  sum  of  the  two  partial  vol¬ 
umes  contained  in  each  vessel.  In  the  specific  example  described  below,  the  rightmost  of 
the  two  cylindrical  vessels  (the  “buffering”  vessel)  has  a  radius  exactly  two  times  that  of 
the  first.  Thus,  any  fluid  added  to  the  leftmost  vessel  (the  “transfer”  vessel)  will  partition 
in  a  ratio  of  1 :4  between  the  transfer  and  buffering  vessels.  The  partial  volume  present  in 
the  transfer  vessel  is  thus  one-fifth  the  total  volume  of  fluid  added  to  the  system.  Like¬ 
wise,  the  partial  volume  present  in  the  buffering  vessel  is  four-fifths  the  total  volume  of 
fluid  added  to  the  system. 

Graphics3D [ { 

Cylinder[{{0,  0,  0},  {0,  0,  5}},  0.5], 

Cylinder[{{3,  0,  0},  {3,  0,  5}},  1], 

Cylinder [{{ 0 . 1 ,  0.1,  0.1},  {3,  0.1,  0.1}},  0.075]}, 

Axes  ^  True,  ImageSize  ^  {300,  300}] 
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□  Formalizing  the  System 

It  is  possible  to  formalize  the  system  below  by  defining  three  functions:  the  transfer  func¬ 
tion  T\{x),  the  buffering  function  ySi(x),  and  the  sigma  function  cri(x). 

The  transfer  function  defines  the  amount  of  fluid  present  in  the  transfer  vessel, 
xi [x_]  :  =  X / 5 

Likewise,  the  buffering  function  defines  the  amount  of  fluid  present  in  the  buffering 
vessel. 


^i[x_]  :  =  4x/5 

The  sigma  function  defines  the  sum  of  the  individual  functions  comprising  the  system. 
ai  [x_]  :  =  xi  [x]  +  lii  [x] 

Once  the  functions  are  defined,  it  is  possible  to  visualize  the  system  by  plotting  the  lines. 

Plot[{ai[x],  /3i[x],  xi[x]},  {x,  0,  100}, 

Axe s Labe 1 ^  {"Total  Volume",  "Partial  Volume"}, 

PlotLegends  "Expressions"] 


-  cr\(x) 

-  ^i(x) 

-  Ti(x) 
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It  is  also  possible  to  visualize  the  system  (in  a  more  intuitive  manner)  by  creating  an  area 
plot  in  which  the  two  partial  volumes  are  indicated  by  the  heights  of  the  transfer  and  buff¬ 
ering  areas  at  any  given  value  of  x. 

Plot[{7;i[x]  ,  ai[x]},  {x,  0,  100}, 

Axe s Labe 1 -»  { "Total\nVolume" ,  "Partial\nVolume" } , 

Filling  ^  Bottom, 

Epilog  ^  { 

Text  [TraditionalForm@HoldForm[^i  [x]  ]  ,  {90,  50}], 

Text  [TraditionalForm@HoldForm[i:i  [x]  ]  ,  {90,  10}], 

Arrowheads [{-. 03,  .03}], 

Arrow[{{80,  0},  {80,  t;i[80]}}], 

Arrow[{{80,  t;i[80]},  {80,  ai[80]}}] 

} 

] 


Partial 

Volume 


□  Defining  the  Buffering  Parameters 

The  definition  of  the  respective  functions  also  provides  a  convenient  means  with  which  to 
quantitate  buffering  action.  This  is  accomplished  in  any  one  of  four  ways,  using  the  three 
available  differentials.  The  buffering  coefficient  measures  the  change  in  the  buffering 
compartment  relative  to  the  total  change. 

bi  =  D[^i[x]]  /D[ai[x]] 


4 

5 
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The  buffering  ratio  Bi  measures  the  change  in  the  buffering  compartment  relative  to  the 
change  in  the  transfer  compartment. 

Bi  =  D[^i [x] ]  /  D[zi [x] ] 

4 

Conversely,  the  transfer  coefficient  ti  measures  the  change  in  the  transfer  compartment  rel¬ 
ative  to  the  total  change. 

ti  =  D[i:i  [x]  ]  /  D[ai  [x]  ] 

1 

5 

The  transfer  ratio  Tj  measures  the  change  in  the  transfer  compartment  relative  to  the  buff¬ 
ering  compartment. 

Ti  =D[7;i[x]]  /D[^i[x]] 

1 

4 

In  this  simple  linear  example — where  the  cross-sectional  area  of  the  buffering  vessel  is 
constant  (i.e.  does  not  vary  with  fluid  level) — the  dimensionless  values  of  total  volumes  b, 
B,  t,  and  T  are  invariant,  equal  to  0.8,  4,  0.20,  and  .0.25,  respectively. 


□  Defining  the  Buffering  Angle 

While  of  great  utility,  the  parameters  B  and  T  are  undefined  in  the  case  of  perfect  buffer¬ 
ing  or  perfect  transfer,  respectively.  To  address  this  issue,  it  is  possible  to  define  an 
alternative  measure,  the  buffering  angle  a,  which  is  capable  of  representing  any  type  of 
buffering  without  the  introduction  of  discontinuities. 
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As  described  in  [1],  the  buffering  angle  a  can  be  defined  by  first  creating  a  three-dimen¬ 
sional  space  curve  of  the  system. 

ParametricPlotSD  [  {D  [x]  ,  D[t:i[x]],  D[^i[x]]},  {x,  0,  100}, 
AxesLabel -»  {x,  y,  z}] 


A  projection  of  this  curve  onto  the  x-y  plane  corresponds  to  the  transfer  function. 

ParametricPlot  [  {D  [x]  ,  D[i:i[x]]},  {x,  0,  100}, 

AxesLabel -»  {x,  y}] 

y 

20  r 

15 
10 
5  - 


A  projection  onto  the  x-z  plane  corresponds  to  the  buffering  function. 

ParametricPlot [ {D [x] ,  D[^i[x]]},  {x,  0,  100}, 
AxesLabel -»  {x,  z}] 


z 
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A  projection  onto  the  y-z  plane,  however,  generates  a  third  curve  that  can  be  used  to  deter¬ 
mine  the  proportions  between  the  individual  changes  in  the  transfer  and  buffering  func¬ 
tions.  This  is  accomplished  by  simply  measuring  the  angle  a\  between  the  y  axis  and  the 
line  formed  by  joining  the  origin  to  any  given  point  on  the  curve. 


Module [{x,  a,  r}  , 

X  =  20; 

a  =  ArcTan[D[i:i  [x]  ]  ,  D[^i[x]]]; 
r  =  Norm[{D[T:i[x]]  ,  D[^i[x]]}]; 

ParametricPlot [ {D [ ti [x] ] ,  D[^i[x]]},  {x,  0,  100}, 
AxesLabel  ^  {y,  z},  ImageSize  ^  80  {1,  3}, 

Epilog  -*  { 

Circle[{0,  0},  r,  {0,  a}], 

Text["ai",  r  {Cos  [1/2  a],  Sin[l/2a]}, 
Background  ^  White] 

} 

] 

] 


z 
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This  calculates  the  buffering  angle  in  degrees  ([1],  Supplement  5). 

Oi  =  If  [ti  ^  1,  1,  -  1]  Nj^ArcCos  1^  Vti^2+bi''2  j]]  /  Degree 

75.9638 

Similar  to  the  buffering  parameters,  is  invariant  in  cases  of  linear  buffering. 

Plot[ai,  {x,  0,  100},  AxesLabel -»  {x,  "ai"}, 

PlotRange  ^  {0,  90}] 


ffi 


80 

60 


40 


20 


20 


40 


60 


80 


100 


It  is  evident  from  the  above  analysis  that  this  methodology  creates  a  system  in  which  per¬ 
fect  buffering  corresponds  to  an  angle  of  90°  and  perfect  transfer  to  an  angle  of  0°.  In  fact, 
any  type  of  buffering  behavior  (moderating/amplifying/inverting/non-inverting)  can  be 
represented  continuously  by  a  finite  angle  between  -45  °  and  135°  ([1],  Supplement  7). 

Importantly,  it  is  also  possible — using  the  same  logic  in  the  reverse  direction — to  use  this 
single  value  to  unambiguously  define  all  of  the  buffering  parameters  describing  the  sys¬ 
tem.  This  can  be  demonstrated  by  defining  a  variable  littlebl  (corresponding  to  the 
buffering  coefficient  b)  that  is  calculated  using  only  the  buffering  angle. 


littlebl = 


Sin[ai °] 

Cos[ai  °]  +Sin[ai  °] 


0 . 8 

This  checks  that  b  and  littleb  are  the  same. 
bi  ::  littlebl 

True 
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Similarly,  the  values  of  the  remaining  three  buffering  parameters  can  be  calculated  using 
only  the  buffering  angle  and  shown  to  be  equal  to  the  corresponding  parameters  calculated 
using  the  relevant  differentials. 

bigBl  =  Sin[ai  °]  /Cos[ai  °] 

4  . 


Bi  ::  bigBl 

True 

littletl  =  Cos [Oi  °  ]  /  (Cos [oi  °  ]  +  Sin [oi  °  ] ) 

0.2 

tl  ::  littletl 

True 

bigTl  =  Cos [oi  ° ]  / Sin [oi  ° ] 

0.25 


Tl  “  bigTl 


True 

Thus,  the  buffering  angle  provides  a  simple,  powerful,  single-value  measure  that  is  capa¬ 
ble  of  capturing  and  communicating  the  underlying  nature  of  the  buffering  relationships 
comprising  the  system. 
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■  A  More  Complex,  Nonlinear  Example 

While  the  simple  linear  scenario  described  above  made  it  possible  to  lay  the  foundation 
for  Schmitt’s  “formal  and  general  approach,”  the  computational  effort  associated  with  its 
application  was  minimal.  Thus,  the  utility  of  applying  these  methods  within  Mathematica 
was  not  fully  exploited.  For  this  reason,  I  next  consider  a  more  complex,  nonlinear  exam¬ 
ple  (i.e.  where  the  cross-sectional  area  of  the  buffering  vessel  is  not  constant).  This  serves 
the  purpose  of  highlighting  the  advantages  of  exploiting  Mathematical  computational 
engine. 


□  Defining  a  Nonlinear  Buffering  System 

Consider  a  nonlinear  system  defined  by  the  functions  listed  below. 

i:2[x_]  :=  0.00395x''11.09746Exp[-2.30471  x] 

02  [x_]  :  =  X 

152  [x_]  :  =  02  [x]  -  1:2  [x] 

As  previously  described,  the  buffering  parameters  can  be  determined  using  the  relevant 
differentials.  However,  in  contrast  to  the  previous  example,  the  computational  burden  asso¬ 
ciated  with  their  calculation  is  no  longer  trivial  (at  least  in  the  absence  of  Mathematica). 

b2  =D[^2[x]]  /D[a2[x]] 

X-0. 00395  0-2- 30471  X  jjll.  0975 


B2  =  D[^2  [x]  ]  /  D[i:2  [x]  ] 

2  53.165  02-30471  X  (x-  0.00395  0-2.30471  x  ^^11.0975) 
jjll.  0975 

t2  =  D[t:2  [x]  ]  /  D[a2  [x]  ] 

0 . 00395  0-2.3O471X  J^IO.  0975 

T2  =  D[i:2  [x]  ]  /  D[^2  [x]  ] 

0 . 00395  0-2.30471  X  jjll.  0975 
X  -  0.00395  0-2.30471  X  jjll.  0975 
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Once  the  functions  are  defined  and  the  system  visualized  through  the  creation  of  the  corre¬ 
sponding  line  and  area  plots,  one  can  clearly  see  that  the  magnitude  of  buffering  action  in 
this  nonlinear  scenario  varies  with  fluid  level  (imagine  a  buffering  vessel  in  the  shape  of 
an  hourglass  as  opposed  to  a  cylinder). 

Plot  [  {1:2  [X]  ,  ^2[X],  02  [x]},  {x,  0,  12}, 

Axe s Labe 1 ^  { "Total \nVolume" ,  "Partial\nVolume" } , 

PlotLegends  -*  "Expressions"] 


Partial 

Volume 


-  Tlix) 

-  Plix) 

-  0-2{x) 


Plot  [{  7:2  [x]  ,  02  [x]},  {x,  0,  12},  Filling  ^  Bottom, 
Axe s Labe 1 -»  { "Total \nVolume" ,  "Partial\nVolume" } , 
Epilog  -»  { 

Text  [HoldForm[/32  [x]  ]  ,  {10,  5}], 

Text  [HoldForm[  1:2  [x]  ]  ,  (5,  1}] 

} 

] 


Partial 

Volume 
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Furthermore,  unlike  the  linear  case,  it  is  also  clear  that  the  buffering  parameters  vary  as  a 
function  of  x. 

Plot[b2,  {x,  0,  12},  PlotRange -»  {0,  2},  AxesLabel  ^  {x,  "b2"}] 


t>2 


Plot[B2,  {x,  0,  12},  PlotRange  ^  {0,  20},  AxesLabel -»  {x,  "B2"}] 


B2 


Plot[t2,  {x,  0,  12},  PlotRange -»  {0,  0.5}, 
AxesLabel -»  {x,  "t2"}] 


ti 
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Plot[T2,  {x,  0,  12},  PlotRange  ^  {0,  1},  AxesLabel  ^  {x,  "T2"}] 


T2 


□  Defining  the  Buffering  Angle  of  a  Nonlinear  System 

As  previously  described,  it  is  also  possible  to  calculate  the  buffering  angle  a2  by  first  creat¬ 
ing  a  three-dimensional  space  curve  of  the  system  and  then  projecting  this  curve  onto  the 
y-z  plane. 

ParametricPlot3D  [  {D  [x]  ,  D[i:2[x]],  D[^2[x]]},  {x,  0,  12}, 
AxesLabel -»  {x,  y,  z}] 
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Module  [{x,  p,  a,  r}  , 

X  =  4; 

p=  {D[r2[x]],  D[^2[x]]}; 
a  =  ArcTan[D[T:2  [x]  ]  ,  D[^2[x]]]; 
r  =  Norm[{D[i:2[x]]  ,  D[^2[x]]}]; 

ParametricPlot  [  {D  [1:2  [x]  ]  ,  D[^2[x]]},  {x,  0,  12}, 
AxesLabel  ^  {y,  z},  ImageSize  ^  100  {1 . 25,  5.5}, 

Epilog  -*  { 

{Orange,  Circle  [{0,  0},  .5r,  (0,  a}]}. 

Text [Row [{"02",  "(",  Style ["X",  Italic],  ")"}], 

.  5  r  (Cos  [1  /  2  a]  ,  Sin[l/2a]},  Background  ^  White]  , 
Thickness [. 02] ,  Orange,  Line[{{0,  0},  p}], 

PointSize [  .  07] ,  Point [p] 

} 

] 

] 


z 
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In  this  instance  the  buffering  angle  a2  varies  as  a  function  of  x,  with  near-perfect  buff¬ 
ering  being  observed  at  both  high  and  low  fluid  levels,  and  minimal  buffering  being  ob¬ 
served  when  X  attains  a  value  of  approximately  4.38. 

02  =  If  [t2  ^  1,  1,  -  1]  N  j^ArcCos  1^  |t2y^  •Vt2''2+b2^2  j]]  /  Degree 

^  ArcCos  [  (O  .  00395  X  2 .11Q2Q-^  ^  j 

j 

0.0000156025  X  2 . 71828-'‘-®°®'‘2  ^  5^20.1949  ^ 
fx  -  0.003  95  X  2.7182  8-2-3°^^^^  5^11.0975^2  ' 


If  [0.00395  e-2-3042ix  5,10.0975  < 


Plot[a2,  {x,  0,  12},  AxesLabel -»  {x,  "02"}^  PlotRange -»  {0,  90}] 


ai 


As  described  previously,  the  buffering  parameters  can  be  determined  using  the  buffering 
angle  and  are  equivalent  to  the  values  calculated  using  the  originally  defined  functions 
(compare  the  plots  below  to  the  plots  for  by,  Bi,  ti,  and  above). 

Iittleb2  =  Sin [02  °  ]  /  (Cos [02  °  ]  +  Sin [02  °  ] )  ; 
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Plot  [littleb2,  {x,  0,  12},  PlotRange  ^  {0,  2}, 
AxesLabel -»  {x,  "littleb2"}] 

littlebl 


bigB2  =  (Sin [02  °]  /  Cos [02  °] ) ; 

Plot[bigB2,  {x,  0,  12},  PlotRange  ^  {0,  20}, 
AxesLabel  ^  {x,  "bigB2"}] 

bigB2 


littlet2  =  Cos [02  ° ]  /  (Cos [02  ° ]  +  Sin [02  °  ] )  ; 
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Plot  [littlet2,  {x,  0,  12},  PlotRange  ^  {0,  0.5}, 
Axe s Labe  1 -»  {x,  "littlet2"}] 

littlet2 


bigT2  =  (Cos [02  °]  /  Sin [02  °] ) ; 


Plot[bigT2,  {x,  0,  12},  PlotRange  ^  {0,  1}, 
AxesLabel ^  {x,  "bigT2"}] 

bigT2 


Thus,  just  as  in  the  case  of  a  linear  system,  plotting  Q'2  over  the  range  of  the  disturbance 
(i.e.  the  addition  of  fluid)  provides  a  simple,  single-value  measure  with  which  to  grasp  the 
inherent  buffering  properties  of  the  system  under  consideration. 
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□  Differentials  versus  Derivatives 

At  this  point  it  is  important  to  note  that,  up  to  now,  only  the  overall  changes  in  the  buffer¬ 
ing  and  transfer  compartments  have  been  considered  (by  comparing  differentials).  How¬ 
ever,  in  the  case  of  nonlinear  systems,  it  is  also  informative  to  consider  the  rate  at  which 
these  values  change  relative  to  one  another  (by  comparing  derivatives).  For  example,  in 
the  same  way  that  it  would  be  important  for  an  investor  to  be  aware  of  the  current  price  of 
gold,  and  how  it  has  changed  in  previous  years,  months,  or  weeks,  it  would  also  be  impor¬ 
tant  for  that  same  investor  to  be  cognizant  of  the  rate  at  which  the  price  is  changing  at  any 
given  instant  in  time  (i.e.  increasing,  decreasing,  or  stable).  This  is  to  say,  the  choice  to 
sell  at  a  given  price  at  a  given  moment  would  be  greatly  influenced  by  knowing  if  the 
price  was  stable  or  on  either  an  upward  or  a  downward  trend.  In  other  words,  it  is  not  only 
of  interest  to  know  how  a  quantity  has  changed,  but  how  that  quantity  is  changing. 

To  determine  the  instantaneous  rates  of  change  in  a  buffered  system,  simply  calculate  and 
plot  the  first  derivatives  of  the  functions  describing  the  system. 

D[^2 ' [X]] 

1-0.043835  e-2. 30471  X  ^10.0975  +  0.0091036  30471  x  ^11.0975 


B[Z2  ' [X]] 

0.04  3835  e-2. 30471  X  ^10.0975  _o. 0091036  ^^11.0975 


D[a2  '  [x]  ] 

1 


Plot  [{ 1:2  '  [x]  ,  ^2  '  [x]  ,  02  '[x]},  {x,  0,  12}, 

Axe s Labe 1 -»  { "Total \nVolume" ,  "Rate  of \nChange"  } , 
PlotLegends  "Expressions"] 


Rate  of 
Change 


Total 

Volume 


-  Ti'ix) 

-  Plix) 

-  0-2  (x) 
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Once  the  functions  are  defined,  the  instantaneous  buffering  parameters  (denoted  by 
adding  an  “i”  to  the  respective  variable  names)  can  also  be  calculated. 

b2i  =  D[^2  '  [X]]  /D[a2  '  [X]] 

1-0.043835  e-2-3047ix  ^10.0975  +  0.0091036  30471  x  ^11.0975 

B2i  =  D[^2  '  [X]  ]  /D[t:2  '  [X]  ] 

(1  -  0.043835  e-2. 30471  X  ^10.0975  +  0.0091036  6^2.30471  x  ^11. 0975 1  ^ 

(0.043835  ^  5^10.0975  _o. 0091036  ^  jjii.0975  | 


t2i  =  D[t:2  '  [x]  ]  /  D[a2  '  [x]  ] 

0.043835  g-2  •  30471  X  ^10.0975  -  0.0091036  g-2  -  30471  x  ^11.0975 

T2i  =  D[T;2'[x]]  /D[^2'[x]] 

(0.04  3835  g-2 -30471  X  jjlO.0975  _  Q  _  QQ9^Q3g  g-2. 30471  x  ^11. 0975  |  ^ 
(1-0.043835  g-2 -30471  X  jjlO.0975  ^  Q  _  QQ9^Q3g  g-2. 30471  x  j^ll.0975  | 

As  is  evident  from  viewing  the  respective  plots  of  these  functions,  this  analysis  allows  one 
to  clearly  identify  the  points  of  maximum  and  minimum  buffering/transfer  rate,  as  well  as 
the  points  at  which  the  buffering/transfer  rates  are  zero. 

Plot[b2i,  {x,  0,  12},  PlotRange  ^  - 1 ,  2}, 

AxesLabel  -»  { "Total \nVolume" ,  HolciForm[^2  '  [x]  /  02  '  [x]  ]  }] 

jSi'W 
0-2' w 
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Plot[t2i,  {x,  0,  12},  PlotRange -»  {-1.5,  1.5}, 

AxesLabel  -»  {  "Total \nVolume" ,  HolciForm[i:2  '  [x]  /  02  '  [x]  ]  }] 

T2'W 


cri'ix) 


Plot[T2i,  {x,  0,  12},  PlotRange^  {-20,  20}, 

AxesLabel  -»  {  "Total \nVolume"  ,  HolciForm[T:2  '  [x]  /  ^2  '  [x]  ]  }] 

Tl'ix) 

p2(x) 


In  this  way,  the  most  critical  junctures  in  the  development  of  the  system  (over  the  range  of 
the  disturbance)  can  be  quickly  and  easily  identified. 
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■  From  the  Abstract  to  the  Practical 

The  analysis  presented  in  the  previous  sections  concludes  my  efforts  to  establish  the  theo¬ 
retical  foundation  of  the  formal  and  general  approach,  as  well  as  its  application  within 
Mathematica.  Next,  in  order  to  illustrate  how  these  abstractions  can  be  directly  applied  to 
the  real  world,  I  consider  an  example  of  non-classical  buffering — namely,  that  of  the 
“genetic  buffering”  of  phenotypic  variation.  In  this  way,  the  ability  to  exactly  and  for¬ 
mally  quantitate  phenomena  that  at  first  glance  might  be  thought  of  as  being  outside  the 
buffering  paradigm  is  practically  demonstrated. 


<=•  Genetic  Buffering 

The  axiomatic  foundation  and  abstract  nature  of  the  “formal  and  general  approach”  allows 
it  to  be  readily  applied  to  any  scenario  in  which  an  arbitrary  quantity  partitions  itself  be¬ 
tween  two  compartments  or  states.  While  suitable  for  the  analysis  of  “classical”  buffering 
phenomena  (the  homeostasis  of  physiological  parameters),  the  approach  is  also  germane 
to  a  variety  of  “non-classical”  buffering  phenomena.  One  such  phenomenon,  which  up  to 
now  has  been  considered  in  only  a  qualitative  manner,  is  that  of  “genetic”  buffering. 
While  many  instances  of  genetic  buffering  have  been  described  and  characterized,  a  well- 
known  and  classic  example  of  the  phenomenon — revolving  around  the  gene  encoding  the 
heat  shock  protein  Hsp90 — will  be  used  to  illustrate  the  real-world  application  of 
Schmitt’s  paradigm. 

■  The  Heat  Shock  Protein  Hsp90 

The  Hsp90  protein  is  an  ATP-dependent  molecular  “chaperone”  that  is  extensively  ex¬ 
pressed  in  organisms  ranging  from  bacteria  to  humans.  It  functions  to  promote  the  proper 
folding  of  a  specific  subset  of  molecular  targets  referred  to  as  its  “client”  proteins.  Interest¬ 
ingly,  many  of  Hsp90’s  client  proteins  are  involved  in  the  process  of  signal  transduction 
and  modulate  developmental  processes  [3].  For  example,  in  the  fruit  fly.  Drosophila 
melanogaster,  mutations  in  the  gene  encoding  Hsp90  result  in  morphological  abnormali¬ 
ties  affecting  the  development  of  the  eye,  legs,  wings,  thorax,  and  bristles  [4]. 

Interestingly,  the  proportion  of  Hsp90  mutant  flies  exhibiting  eye  defects  (in  a  line  prone 
to  such  abnormalities)  increases  dramatically  as  the  temperature  rises.  In  contrast,  the  pro¬ 
portion  of  control  flies  (expressing  the  normal  or  “wild-type”  version  of  the  Hsp90  pro¬ 
tein)  also  increases  in  response  to  increased  temperature,  but  not  nearly  to  the  same  degree 
as  in  the  mutant  flies  [4].  The  gene  encoding  the  wild-type  version  of  the  Hsp90  protein  is 
thus  said  to  “buffer”  the  appearance  of  the  defective  eye  trait. 
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■  Formally  Quantitating  Hsp90-Medlated  Buffering  of  Phenotypic  Variation 

It  is  possible  to  analyze  these  experimental  results  in  more  detail  by  importing  the  raw 
data  (see  Figure  8  in  [4])  from  the  original  publication  in  the  form  of  two  arrays.  These 
two  arrays,  corresponding  to  the  wild-type  and  mutant  fly  lines,  relate  temperature  to  the 
penetrance  of  the  mutant  phenotype  (i.e.  the  percentage  of  flies  in  the  line  displaying  the 
eye  defect). 

wildtypeHsp90  =  {{18,  0},  {21,  0},  {25,  0},  {26,  1}, 

{27,  3},  {30,  19},  {32,  16}} 

{{18,  0),  {21,  0},  {25,  0), 

{26,  1},  {27,  3},  {30,  19},  {32,  16}} 


mutantHsp90  =  {{18,  2},  {21,  10},  {25,  10},  {26,  25}, 
{27,  45},  {30,  70},  {32,  82}} 

{{18,  2},  {21,  10},  {25,  10}, 

{26,  25},  {27,  45},  {30,  70},  {32,  82}} 


Plotting  the  data  illustrates  that  the  mutant  line  is  indeed  more  sensitive  to  increasing  tem¬ 
perature  than  the  wild-type  line. 

Show[ListPlot [mutantHsp90,  PlotStyle  Red, 

PlotMarkers  ^Automatic, 

Axe s Labe 1  ^  {"Temperature  (°C)",  "Affected  (%)"}, 

PlotRange  -»{-l,  100}], 

ListPlot [wildtypeHsp90,  PlotStyle  -* Darker@Green, 

PlotMarkers  Automatic]  ] 


Affected  (%) 

too  r 


80 

60 


40 


20 


•  • 

•  • 

'  e  — . 9 ' 9  '  ^ . — ^  Temperature  {°C) 

20  22  24  26  28  30  32 
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To  proceed  with  the  application  of  the  formal  and  general  approach,  we  now  create  two 
curves  that  fit  the  data  using  the  FindFit  function. 

FindFit  [wildtypeHsp90,  jExp[-kX],  {j,  k},  X, 

Method  NMinimize] 

{  j  ^  0 . 0031044,  k^-  0. 272652} 


FindFit  [mutantHsp90,  jExp[-kX],  {j,  k},  X, 

Method  NMinimize] 

{j^O. 193965,  k^-  0. 191242} 

With  the  above  constants,  the  curves  describing  the  behavior  of  the  wild-type  and  mutant 
lines  can  be  generated. 

wildtypecurve [x_]  := 

jExp[-kx]  /.  FindFit  [wildtypeHsp90,  jExp[-kX], 

{j,  k}  ,  X,  Method  ^  NMinimize] 

mutantcurve [x_]  := 

jExp[-kx]  /.  FindFit  [mutantHsp90,  jExp[-kX],  {j,  k}, 

X,  Method  NMinimize] 


wildtypecurve [x] 

0 . 0031044  e°-2^2652  x 


mutantcurve [x] 

0 .193965  e°-i®^242x 
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The  appropriateness  of  the  fit  can  then  be  inspected  by  plotting  the  curves  together  with 
the  raw  data. 

Show[ 

ListPlot [mutantHsp90,  PlotStyle  ^  Red, 

PlotMarkers  Automatic, 

AxesLabel  ^  {"Temperature  (°C)",  "Affected  (%)"}], 

ListPlot  [wildtypeHsp90,  PlotStyle  ^  Darker@Green, 

PlotMarkers  Automatic]  , 

Plot [Evaluate [ {mutantcurve [x] ,  wildtypecurve [x] } ] , 

(x,  18,  32},  PlotStyle  ^  {Red,  Darker@Green} ] 


With  these  results  in  hand,  it  is  now  possible  to  define  both  the  buffering  and  transfer  func¬ 
tions  and  then  to  visualize  the  system  using  an  area  plot. 

^3 [x_]  :  =  Evaluate [wildtypecurve [x] ] 


/^3  [^_]  •  =  Evaluate  [mutantcurve  [x]  ] 


03  [X_]  :  =  Z3  [x]  +  ^3  [x] 
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Plot  [  {ts  [x]  ,  03  [x]},  {x,  18,  32},  Filling  ^  Bottom, 
AxesLabel ^  { "Temperature  ( °C) " ,  "Affected  (%)"}, 
Epilog  ^  { 

Text [HoldForm[^3  [x]  ]  ,  {31,  40}], 

Text  [HoldForm[  1:3  [x]  ]  ,  {31,  7}] 

} 

] 


Affected  (%) 


Temperature  (°C> 


As  is  evident  from  inspecting  this  plot,  an  increase  in  temperature  translates  into  only  a 
small  increase  in  the  proportion  of  flies  with  eye  defects  in  the  wild-type  line.  However,  in 
the  mutant  line,  a  much  more  pronounced  increase  is  observed.  We  can  thus  infer  that  the 
expression  of  the  wild-type  Hsp90  protein  buffers  the  effect  of  temperature  on  the  pene¬ 
trance  of  the  eye  defect.  In  other  words,  since  the  “transfer”  of  the  eye  defect  phenotype  is 
greater  in  the  mutant  line,  the  wild-type  Hsp90  protein  must  possess  the  abstract  capacity 
of  diverting  or  “soaking  up”  this  disturbance  (into  an  equally  abstract  buffering  compart¬ 
ment)  so  that  the  “transfer”  of  the  eye  defect  is  diminished.  Importantly,  using  the  formal 
and  general  approach,  we  can  now  provide  an  exact  quantitative  measure  of  the  buffering 
capacity  of  the  wild-type  allele  by  calculating  and  plotting  the  four  buffering  parameters 
together  with  the  buffering  angle  (below). 

b3  =  D[/33[x]]  /D[a3[x]] 

0 .193965 

0.193965  +  0.0031044  e°-2’^2652  x 
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Plot[b3,  {x,  18,  32},  PlotRange  ^  .  8,  1}, 

AxesLabel ^  {"Temperature  (°C)",  "bs"}] 


Temperature  (°C) 


B3  =  D[^3[x]]  /D[t:3[x]] 

62.4808 


Plot[B3,  (x,  18,  32},  PlotRange -*  {4,  16}, 
AxesLabel -»{ "Temperature  (°C)",  "B3"}] 


t3  =  D[i:3  [x]  ]  /  D[a3  [x]  ] 

0 .0031044  e“'2’^2652  x 

0.193965  e“- 151242  X  +  0.0031044  e°  •  212652  x 
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T3  =  D[t:3  [x]  ]  /  D[^3  [x]  ] 

0 .0160049 

Plot[T3,  {x,  18,  32},  AxesLabel ^  {"Temperature  (°C)",  "T3"}] 
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As  shown  above,  the  buffering  coefficient  can  be  calculated  to  range  from  -0.935  at  18 
°C  to  -0.822  at  32  °C.  While  these  calculations  are  now  trivial,  the  critical  biological  ques¬ 
tion  of  how  to  interpret  these  calculations  is  more  challenging.  As  discussed  by  Schmitt 
([1],  Supplements  6  and  7),  one  way  to  proceed  is  to  consider  the  conceptual  overlap  be¬ 
tween  the  formal  and  general  approach  and  probability  theory.  In  essence,  the  buffering 
coefficient  is  analogous  to  a  probability.  In  the  same  way  that  a  probability  measures 
the  proportion  of  a  part  to  a  whole  (i.e.  the  number  of  successful  events  to  the  total  num¬ 
ber  of  events),  so  does  the  buffering  coefficient  (by  measuring  the  fractional  change  in  the 
buffering  compartment  relative  to  the  whole  compartment).  The  critical  realization  is  that 
both  quantities  are  measured  on  a  relative  scale  that  is  normalized  to  1.  Interpreting  the 
graph  in  this  way,  one  can  say  that  the  “probability”  of  buffering  decreases  from  -0.935  to 
-0.822  from  18  °C  to  32  °C. 

Using  similar  logic,  one  can  say  that,  in  essence,  the  buffering  ratio  is  analogous  to  the 
“odds”  of  buffering.  In  the  same  way  that  odds  measure  the  proportion  between  two  parts 
of  a  whole  (by  comparing  events  to  non-events),  so  does  (by  comparing  the  fractional 
change  in  the  buffering  compartment  to  the  fractional  change  in  the  transfer  compart¬ 
ment).  Again  the  critical  realization  is  that  both  quantities  are  measured  on  a  scale  with 
equal  intervals  and  an  absolute  zero.  Interpreting  the  data  in  this  way,  one  can  say  that  the 
“odds”  of  buffering  decrease  from  -14.4  to  -4.6  from  18  °C  to  32  °C.  Similarly,  one  can 
consider  and  T-^  to  be  analogous  to  the  “probability”  of  transfer  and  the  “odds”  of  trans¬ 
fer,  respectively. 

Lastly,  the  buffering  angle  provides  information  regarding  not  only  the  magnitude  of 
buffering  but  also  the  class  of  buffering  behavior.  As  discussed  by  Schmitt  ([1],  Supple¬ 
ment  7),  angles  between  90°  and  0°  correspond  to  non-inverting  moderation,  with  angles 
of  90°  corresponding  to  perfect  buffering  and  0°  to  zero  buffering.  Thus,  the  simple  inspec¬ 
tion  of  the  plot  of  a-^  reveals  that  the  buffering  capacity  of  wild-type  Hsp90  is  very  strong 
(near  perfect  at  -86°)  but  decreases  slightly  with  increasing  temperature  (to  a  value  of 
-78°).  In  any  event — and  however  the  data  is  interpreted — the  analysis  presented  above 
clearly  demonstrates  the  ease  with  which  the  buffering  capacity  of  a  genetic  determinant 
can  be  formally  and  exactly  quantified. 
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■  Conclusion 

A  formal  conceptualization  of  genetic  buffering  was  first  put  forward  by  famed  geneticist 
and  Nobel  laureate  Leland  Hartwell  [5,  6].  In  these  publications,  Hartwell  highlights  the 
importance  of  buffering  relationships  in  determining  phenotype  in  outbred  organisms. 
Subsequent  research  both  in  yeast  and  in  the  roundworm,  Caenorhabditis  elegans,  has 
strongly  supported  his  assumptions  [7,  8].  Thus,  in  addition  to  increasing  our  basic  under¬ 
standing  of  genetic  networks,  these  results  also  raised  the  theoretical  question  of  whether 
or  not  it  might  be  possible  to  provide  an  exact  quantitative  measure  of  a  genetic  determi¬ 
nant’s  buffering  activity. 

As  described  for  the  first  time  above,  Schmitt’s  “formal  and  general  approach”  provides  a 
powerful  means  with  which  to  accomplish  this  goal.  Through  the  calculation  of  the  four 
buffering  parameters  b,  B,  t,  and  T  and  the  buffering  angle  a,  it  is  indeed  possible  to 
provide  an  exact  quantitative  measure  of  the  buffering  activity  or  “power”  of  any  genetic 
determinant  in  response  to  an  environmental  disturbance.  Furthermore,  the  activity  can  be 
measured  using  either  a  relative  scale  normalized  to  one  {b  and  t)  or  an  absolute  scale  with 
equal  intervals  (B  and  T),  or  distilled  to  a  single  value  measure  a,  capable  of  describing 
the  system  in  its  entirety. 

In  conclusion,  it  must  be  noted  that  the  application  of  the  formal  and  general  approach  be¬ 
comes  practical  only  when  used  in  conjunction  with  a  sophisticated  computational  engine 
that  is  capable  of  readily  analyzing  the  relevant  data.  As  demonstrated  in  this  article,  Math- 
ematica  provides  such  an  engine.  The  application  of  the  formal  and  general  approach 
using  the  Wolfram  Language  thus  provides  a  powerful  tool  with  which  to  quantitatively 
analyze  both  classical  and  non-classical  buffering  phenomena,  irrespective  of  discipline. 


■  References 


[1]  B.  M.  Schmitt,  'The  Quantitation  of  Buffering  Action  i.  A  Formai  &  Generai  Approach,”  Theo¬ 
retical  Biology  and  Medical  Modelling,  2(8),  2005.  doi:1 0.1 1 86/1742-4682-2-8. 

[2]  B.  M.  Schmitt,  ‘The  Quantitation  of  Buffering  Action  ii.  Appiications  of  the  Formai  &  Generai 
Approach,”  Theoretical  Biology  and  Medical  Modelling,  2(9),  2005. 

doi:1 0.1 186/1742-4682-2-9. 

[3]  M.  Taipaie,  D.  F.  Jarosz,  and  S.  Lindquist,  “FiSP90  at  the  Fiub  of  Protein  Fiomeostasis: 
Emerging  Mechanistic  insights,”  Nature  Reviews  Molecular  Cell  Biology,  11(7),  2010 
pp.  515-528.  doi:10.1038/nrm2918. 

[4]  S.  L.  Rutherford  and  S.  Lindquist,  “F1sp90  as  a  Capacitor  for  Morphoiogicai  Evoiution,”  Na¬ 
ture,  396,  1998  pp.  336-342.  doi:1 0.1 038/24550. 

[5]  J.  L.  Fiartman,  B.  Garvik,  and  L.  Fiartweii,  “Principies  for  the  Buffering  of  Genetic  Variation,” 
Science,  291(5506),  2001  pp.  1001-1004. 
www.sciencemag.org/content/291/5506/1001  .abstract. 

[6]  L.  Fiartweii,  “Robust  interactions,”  Science,  303(5659),  2004  pp.  774-775. 
doi:1 0.1 126/science.1 094731. 


The  Mathematica  Journal  17  ©2015  Wolfram  Media,  Inc. 


30 


Jim  Karagiannis 


[7]  A.  H.  Tong,  G.  Lesage,  G.  D.  Bader,  H.  Ding,  H.  Xu,  X.  Xin,  J.  Young,  G.  F.  Berriz,  R.  L. 
Brost,  M.  Chang,  Y.  Chen,  X.  Cheng,  G.  Chua,  H.  Friesen,  D.  S.  Goldberg,  J.  Haynes,  C. 
Humphries,  G.  He,  S.  Hussein,  L.  Ke,  N.  Krogan,  Z.  Li,  J.  N.  Levinson,  H.  Lu,  P.  Menard,  C. 
Munyana,  A.  B.  Parsons,  O.  Ryan,  R.  Tonikian,  T.  Roberts,  A.  M.  Sdicu,  J.  Shapiro,  B. 
Sheikh,  B.  Suter,  S.  L.  Wong,  L.  V.  Zhang,  H.  Zhu,  C.  G.  Burd,  S.  Munro,  C.  Sander,  J.  Rine, 
J.  Greenblatt,  M.  Peter,  A.  Bretscher,  G.  Bell,  F.  P.  Roth,  G.  W.  Brown,  B.  Andrews,  H. 
Bussey,  and  C.  Boone,  “Global  Mapping  of  the  Yeast  Genetic  Interaction  Network,’’  Science, 
303(5659),  2004  pp.  808-813.  doi:10.1 126/science.1091317. 

[8]  B.  Lehner,  C.  Crombie,  J.  Tischler,  A.  Fortunate,  and  A.  G.  Fraser,  “  Systematic  Mapping  of 
Genetic  Interactions  in  Caenorhabditis  Elegans  Identifies  Common  Modifiers  of  Diverse  Sig¬ 
naling  Pathways,’’  Nature  Genetics,  38(8),  2006  pp.  896-903.  doi:10.1038/ng1844. 

J.  Karagiannis,  “The  Quantitation  of  Non-classical  Buffering,”  The  Mathematica  Journal,  2015. 
dx.doi.org/doi:10.3888/tmj.17-3. 


About  the  Author 

Dr.  Jim  Karagiannis  is  an  associate  professor  at  the  University  of  Western  Ontario  in  Lon¬ 
don,  Ontario,  Canada.  His  research  makes  use  of  the  model  eukaryote  Schizosaccha- 
romyces  pombe  and  explores  the  complex  post-translational  modifications  that  take  place 
on  the  carboxy-terminal  domain  (CTD)  of  the  largest  subunit  of  RNA  polymerase  II. 
Through  an  empirical  examination  of  the  informational  properties  and  regulatory  potential 
of  the  CTD,  Dr.  Karagiannis  hopes  to  decipher  the  “programming  language”  used  by 
eukaryotes  to  control  aspects  of  gene  expression. 

Jim  Karagiannis 

Department  of  Biology 

The  University  of  Western  Ontario 

London,  ON 

Canada,  N6A-5B7 

jkaragia  @  uwo.  ca 


The  Mathematica  Journal  17  ©201 5  Wolfram  Media,  Inc. 


